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Abstract 

The equilibrium phase behavior of a binary mixture of charged colloids and neutral, non- 
adsorbing polymers is studied within free-volume theory. A model mixture of charged hard-sphere 
macroions and ideal, coarse-grained, effective-sphere polymers is mapped first onto a binary hard- 
sphere mixture with non-additive diameters and then onto an effective Asakura-Oosawa model 
[S. Asakura and F. Oosawa, J. Chem. Phys. 22, 1255 (1954)]. The effective model is defined by a 
single dimensionless parameter - the ratio of the polymer diameter to the effective colloid diameter. 
For high salt-to-counterion concentration ratios, a free- volume approximation for the free energy is 
used to compute the fluid phase diagram, which describes demixing into colloid-rich (liquid) and 
colloid-poor (vapor) phases. Increasing the range of electrostatic interactions shifts the demixing 
binodal toward higher polymer concentration, stabilizing the mixture. The enhanced stability is at- 
tributed to a weakening of polymer depletion-induced attraction between electrostatically repelling 
macroions. Comparison with predictions of density-functional theory reveals a corresponding in- 
crease in the liquid-vapor interfacial tension. The predicted trends in phase stability are consistent 
with observed behavior of protein-polysaccharide mixtures in food colloids. 
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I. INTRODUCTION 



Mixtures of colloidal particles and free (non-adsorbing) polymer coils dispersed in a solvent 
are among the most intensively studied soft matter system o 1 ! 2 ! 3 ! 4 ! 5 ! 6 . The conceptual analogy 
between colloids and atoms, similarities in thermodynamic phase behavior between colloidal 
suspensions and atomic systems, and the relative ease of tuning polymer-induced effective 
colloidal interactions, make colloid-polymer mixtures valuable model systems for probing 
connections between microscopic interparticle interactions and macroscopic properties in a 
variety of materials. The observed phases, distinguished by composition and structural order, 
include most of the equilibrium phases familiar in simple molecular systems, e.g., vapor, 
liquid, crystal, as well as nonequilibrium states, such as glasses and gels^. Complementing 
their fundamental importance, colloid-polymer mixtures have diverse industrial applications, 
e.g., to coatings, petroleum products, pharmaceuticals, and many foods, where polymer 
additives are used to control phase stability and rheological properties. 

Addition of free polymer can substantially modify effective interactions between colloidal 
particles through the mechanism of entropic depletion. When two colloidal surfaces approach 
to a separation closer than the typical diameter of a polymer coil, the entropic cost to the 
coil of distorting its average spherical conformation tends to exclude the polymer. The 
resulting depletion of polymer from the space between the colloids creates an imbalance in 
polymer osmotic pressure that can induce effective attractions between colloids. Attractions 
of sufficient range and strength, depending on the relative size and concentration of polymer, 
can drive bulk demixing into colloid-rich (liquid) and colloid-poor (vapor) phases. 

The first, and conceptually simplest, statistical mechanical model to qualitatively de- 
scribe the polymer depletion mechanism and the associated phenomenon of depletion-driven 
demixing was the model of Asakura and Oosawa£, developed independently by Vrij*- The 
Asakura-Oosawa (AO) model regards the colloidal particles as hard spheres and the polymer 
coils as effective ideal spheres - mutually noninteracting, but interacting with the colloids 
as hard spheres. Neglect of polymer-polymer interactions can be reasonably justified for 
polymers in a theta solvent. The phase diagram corresponding to the AO model has been 
determined by a variety of methods, including thermodynamic perturbation theory^, free- 
volume theor y 11 ! 12 , density-functional theory^ 3 -, and Monte Carlo simulatio n 14 ! 15 ! 16 . Recent 
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effort has also been devoted to calculating the liquid-vapor interfacial tension of phase- 
separated colloid-polymer mixture a 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 i 26 . 

Although certain topologically constrained polymers (e.g., stars and microgels) are close 
to spherical in shape, linear-chain polymers are random walks whose shapes are spherical 
only on average. To explore the significance of nonspherical conformations, more explicit 
segmented-chain polymer models have been studied via Monte Carlo simulatio n 27 i 28 i 29 and 
integral-equation theory 3 ^ Effective polymer-polymer interactions have been modeled by 
combining Monte Carlo and integral-equation methods within a "polymers as soft colloids" 
framework^ 9 -. In recent work, we examined fluid-fluid demixing within several variations of 
the classic AO model, incorporating a third component (hard needles)^, polymer-solvent 
interactions* 3 - 2 -, polymer-polymer interactions** 3 -, and colloid- induced polymer compression 34 . 
The latter study incorporated intrinsic polydispersity in the polymer radius of gyration. 
Influences on phase behavior of colloid polydispersity 3 ^ and of polymer chain length poly- 
dispersity 3 ^ 7 - have also been studied theoretically. 

While neutral colloid-polymer mixtures have been widely investigated, much less atten- 
tion has been devoted to mixtures of charged species. In the case of charged colloids or 
charged polymers (polyelectrolytes) , electrostatic repulsions between macroions, screened 
by counterions and salt ions in solution, compete with depletion-induced attractions and 
can modify phase behavior. In pioneering experiments** 8 -!** 9 -, phase separation was observed 
in mixtures of charged colloids and neutral polymers. Gast et a/.— interpreted their ob- 
servations by applying thermodynamic perturbation theory to an effective one-component 
model, incorporating electrostatic repulsion and depletion-induced attraction into an effec- 
tive pair potential between colloids. Subsequent measurements of force profiles in mixtures 
of charge-stabilized colloidal oil-in-water emulsion droplets and ionic surfactant micelles 4 ^ 
and in mixtures of colloidal particles and charged macromolecules 4 ^ directly demonstrated 
the potential for electrostatic interactions to modify depletion forces. More recently, liquid- 
vapor separation and gelation were observed in mixtures of charged colloids and ionic worm- 
like micelles* 4 ^ and depletion potentials induced by charged rods were measured in colloidal 
rod-sphere mixtures 4 * 3 -. In a recent theoretical study, liquid-state integral-equation methods 
were used to model the structure and phase behavior of mixtures of charged colloids and 
polyelectrolytes 44 . 

3 



The purpose of the present paper is to propose, for a simple model of a charged-colloid 

- neutral-polymer mixture, an alternative theoretical approach that, in contrast to Ref. |3j| 
treats the two components on an equal footing. The theory is based on mapping the binary 
mixture onto an effective AO model, governed only by excluded volume interactions, and 
applying free- volume theory^ and classical density-functional theory^. Within this concep- 
tual approach, we explore the qualitative influence of colloidal charge on demixing and find 
that increasing the range of electrostatic interactions can significantly stabilize the mixture 
against demixing and correspondingly increase the liquid-vapor interfacial tension. 

The remainder of the paper is organized as follows. The model system is defined in 
Sec. IHl The mapping onto the effective AO model and the free-volume theory are described 
in Sec. 11111 Results for the fluid-fluid demixing phase diagram and interfacial tension are 
presented and discussed in Sec. II VI Finally, Sec. |V| closes with a summary and conclusions. 

II. MODEL 

The system of interest comprises charged colloidal particles, their dissociated counterions, 
and free polymer coils, all dispersed in an electrolyte solvent. The multi-component mix- 
ture of macroions, microions (counterions and salt ions), polymers, and solvent molecules 
spans a range of length and time scales, presenting severe challenges to explicit modeling 
approaches. To reduce the system to a tractable model, we make several simplifying, yet 
realistic, assumptions. (1) The colloids are modeled as monodisperse charged hard spheres 

- a close approximation to many real synthetic suspensions. (2) The microions are modeled 
as monovalent point charges and subsumed into effective electrostatic interactions between 
macroions, as described below. (3) The polymer coils are represented, as in the coarse- 
grained AO model, as effective spheres. Nonspherical conformations, disfavored by lower 
conformational entropy, can be reasonably neglected when the polymer coils are comparable 
or smaller in size than the colloids (colloid limit). (4) Finally, the solvent is treated as a di- 
electric continuum, characterized by a dielectric constant (primitive model of electrolytes^). 
An experimental system that closely resembles our model system would be charge-stabilized 
synthetic polystyrene or silica microspheres dispersed in an aqueous electrolyte together with 
a non-adsorbing, nonionic, water-soluble polymer, such as hydroxyethylcellulose (HEC) 38 ' 39 
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or polyethylene oxide (PEO). 

Interactions between colloidal particles include steric and electrostatic repulsions and van 
der Waals attractions. Here we assume hard-sphere steric interactions and ignore van der 
Waals interactions, which is valid for particles index- matched to the solvent. Electrostatic 
interactions between colloidal macroions result from bare Coulomb repulsion and screening 
by surrounding micro ions. The classic theory of Derjaguin, Landau, Verwey, and Overbeek 
(DLVO)^ predicts that, in a dilute suspension, two colloidal macroions at center-to-center 
separation r interact via a pair potential of screened-Coulomb (Yukawa) form: 



v cc (r) 



r < 2R r 



1 + kR c j r ' r ^ 2R o 



Z 2 e 2 / 6 kR c \ 2 e -KT 



where Z is the effective macroion valence, e the proton charge, e the solvent dielectric 
constant, k the Debye screening constant (inverse screening length), and R c the colloid 
radius. 

More recent statistical mechanical approaches proceed by formally integrating out from 
the partition function the microion degrees of freedom to map the macroion-microion mix- 
ture onto an effective one-component system governed by effective interaction a 47 i 48 i 49 i 50 i 51 i 52 , 
determined by the distribution of microions around the macroions. The microion distri- 
bution depends on the response - in general nonlinear - of the microions to the macroion 
charge density. For sufficiently weakly charged macroions, the microion response can be 
approximated as linear. A further neglect of microion correlations, which is reasonable for 
monovalent counterions, then recovers the DLVO form of effective pair potential [Eq. ((TJ] 
with a density- dependent screening constant of the form 



K= i — — ' (2) 

where p c and p s are, respectively, the number densities of colloidal macroions and salt ion 
pairs, k-Q is Boltzmann's constant and T is the absolute temperature. The nonzero size of 
the microions may modify the microion distribution and thus screening effect, but should 
not qualitatively change the general repulsive form of the effective pair potential assumed 
here. In practice, since most water-soluble polymers are less polarizable than water, e may 
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be best interpreted as an effective dielectric constant of the polymer-solvent mixture and k 
as an effective screening constant. For the purposes of this paper, the essential point is that 
the range of electrostatic repulsion, governed by k, can be widely tuned by adjusting the 
ionic strength (salt concentration) of the electrolyte solvent. 

As a natural byproduct of the one-component mapping, the total free energy contains, in 
addition to a pair-interaction contribution, a one-body volume energ y 47 i 48 i 49 i 50 i 51 i 52 E, given 
by 

^ = (Z A+ft ) {>n[(Zp c + ,,)A3] - 1}+,, {1 „(M») - 

(3) 

where V is the total volume of the system, A is the thermal wavelength of the microions and 
Ab = e 2 /(ekftT) is the Bjerrum length. In a straightforward physical interpretation, the first 
two terms in Eq. (JHJ) account for the microion entropy, the third term is the interaction of a 
macroion with its own cloud of counterions, and the final term results from charge neutrality. 
Because the volume energy depends nontrivially on colloid density, it must, in general, be 
included in the free energy. At relatively high salt concentrations (p s 3> Zp c ), however, the 
counterions are dominated by salt ions and the volume energy can be neglected. We restrict 
our considerations to this parameter regime. 

Within the coarse-grained spherical polymer model, colloid-polymer interactions are de- 
scribed by a simple excluded-volume pair potential: 

I oo, r < R c + R p 
Vc P (r) = { (4) 
{ 0, r > R c + Rp, 

where r is now the macroion-polymer center-to-center distance and R p is the polymer ra- 
dius of gyration. For simplicity, we assume the solvent to be near its theta temperature 
for the polymer—, allowing the polymer to be reasonably modeled as ideal (mutually non- 
interacting), i.e., v pp (r) = 0, for all r. In practice, the strength of polymer-polymer interac- 
tions will depend on the properties of specific polymers and solvents. 

Within the above assumptions, the model system is now completely characterized by the 
colloid-colloid interaction [Eqs. (JTJ) and (J2J)] and the polymer-to-colloid size ratio q = R p /R c . 
The thermodynamic states of the system are specified by the bulk number densities pi, or 
equivalently the volume fractions, rji = (47r/3)pji??, of species % — c,p. 
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III. THEORY 

A. Mapping onto the Asakura-Oosawa Model 

To explore thermodynamic properties of the model colloid-polymer mixture, we seek an 
approximation for the free energy of the system. We proceed by first constructing a mapping 
onto a simpler model system. A similar approach has been developed recently by Tuinier— . 
Assuming the screened- Coulomb effective pair potential between colloids [Eq. (0)] to be 
relatively steeply repulsive (kR c > 1), the corresponding contribution to the free energy 
may be mapped, with reasonable accuracy, onto the free energy of an effective hard-sphere 
syste m 45 ! 56 interacting via an effective hard-sphere pair potential, 



where R' c is the effective hard-sphere colloid radius. A reasonable first estimate of the effective 
radius is obtained from a simple thermal criterion, 



according to which colloids tend not to approach closer than a distance at which their 
interaction energy is comparable to the typical thermal energy. 

Assuming high salt-to-counterion concentration ratios (p s ^> Zp c ), the screening constant 
k, and thus v cc (r) and R' c , are practically independent of density. In this same limit, the 
volume energy [Eq. ©] contributes to the free energy density a term that is only linear 
in colloid density and therefore irrelevant for phase behavior. To demonstrate that in this 
high-salt-concentration regime the effective colloid radius may still significantly exceed the 
bare radius, consider the typical example of particles of bare radius R c = 50 nm and effective 
valence Z = 500 suspended at volume fraction r] c = 0.01 in a 1 mM aqueous salt solution 
at room temperature. In this case, p s ^> Zp c = 0.016 mM, yet the effective colloid radius is 
estimated [from Eq. (jUJ)] to be R' c ~ 67 nm, i.e., about 30% greater than the bare radius, a 
difference that also well exceeds typical colloid polydispersities. 

The original system has thus far been mapped onto an effective binary hard-sphere mix- 
ture, governed by the pair interactions of Eqs. (@J) and (j3J). As a measure of the relative range 
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of electrostatic and excluded-volume repulsions, it is convenient to define the electrostatic 
colloid size ratio, £ = R' c /R c > 1, as the ratio of the effective and bare colloid radii. The 
model system is now competely specified by two dimensionless parameters, £ and q. Since 
the polymer species is ideal, this model resembles the AO model, with the colloid diameter 
merely rescaled. However, because colloids and polymers repel at a range of (R c + R p ), 
rather than (R! c + Rp), the hard-sphere diameters are no longer additive and the rescaling 
is nontrivial. Nevertheless, as Fig. ^ illustrates, this non-additive mixture can be further 
mapped exactly onto an effective additive AO model by introducing a fictitious effective 
polymer radius R' p , defined via 

R' c + R' p = R c + R P (7) 

so as to yield the correct colloid-polymer excluded-volume interaction [Eq. (flj)]. 

The final effective AO model is characterized by a single dimensionless parameter, namely 
the effective polymer-to-colloid size ratio 

q' = R' p /R' c = (q-£ + l)/Z, for £ < 1 + q. (8) 

As shown in Fig. |21 the rescaled polymer-to-colloid size ratio is always smaller than the 
true ratio (</ < q) and with increasing electrostatic colloid size ratio, q' decreases, i.e., the 
polymers grow effectively smaller relative to the colloids. For £ > 1 + q, there is no polymer 
depletion, and so then q' = 0. The thermodynamic states are in turn specified by effective 
volume fractions 

rf = yPX 3 = fVc (9) 

and 

Air 

v' P = -jP P R' P 3 = [!-(£- (10) 

where we have used the relation R' p / R p — 1 — (£ — l)/q < 1, obtained from Eq. (J2J). 
B. Free- Volume Theory 

Having mapped the original mixture of charged colloids and neutral polymers onto an 
effective AO model, we now consider the Helmholtz free energy, from which all equilibrium 
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thermodynamic properties may be determined. The total free energy, F = + F ex , sep- 
arates conveniently into an ideal-gas term F^, which is independent of interactions, and 
an excess term F ex , which depends entirely on interactions. For a homogeneous fluid, the 
ideal-gas free energy density is given exactly by 

PF d /V = p c [ln(p c A c 3 ) - 1] + p P [HP P K) ~ 1], (11) 

where (3 = l/k^T and A c and A p are the thermal wavelengths of the colloid and polymer 
species. For the excess free energy, we apply the mean-field free-volume theory of Lekkerk- 
erker et a/.—, which predicts fluid-fluid phase separation for the original AO model in good 
agreement with simulatio n 14 ' 15 ' 23 ' 24 . Within this approach, the fluid excess free energy den- 
sity is approximated by 

(3FJV = Pfastfc) ~ P P lna( V ' c ,q>), (12) 

where 4>hs(v' c ) * s the excess free energy density of the effective one-component hard-sphere 
fluid and cv(r)' c ,q ; ) is the polymer free-volume fraction, i.e., the fraction of the total volume 
accessible to the polymer centers (not excluded by the colloids). The free- volume fraction 
is related to the excess chemical potential of the polymer via Widom's particle insertion 
method 57 : 

A*p,ex = —k B T\n(exp(—AU/kBT)) = — /c B Tlna, (13) 

where the angular brackets denote an ensemble average over colloid configurations and poly- 
mer positions, and AU is the change in potential energy upon insertion of a polymer - 
here infinite if the polymer overlaps a colloid and zero otherwise. An approximation for 
a is obtained by noting that /i p>ex also equals the reversible work required to produce a 
spherical cavity, of radius R' p , in a fluid of hard sheres of radius R' c . The scaled-particle 
approximatio n 58 ' 59 for /x P)ex then yields 

a(r}' c ,q') = (1 - rf c ) exp(-Ar - Bx 2 - Cx 3 ), (14) 

with x = rj'J(l - r}' c ), A = q ri + 3q' 2 + 3q', B = 3q' 3 + 9q' 2 /2, and C = 3q' 3 . Within the same 
approximation, the excess free energy of the colloids is given by 

P***™ ~ 8vr^3(l-^) 2 • ( } 



The pressure resulting from Eq. (|T5|) coincides with the compressibility equation of state 
following from the exact solution of the Percus-Yevick integral equation for hard spheres 60 ' 61 . 
Taken together, Eqs. (|11|) - (|15|) provide an approximation for the fluid free energy, from which 
the fluid-fluid demixing phase diagram can be computed. 

IV. RESULTS AND DISCUSSION 

Based on the approximate free energy of Eqs. I)11 J) -(|15 )) . we have performed a coexistence 
analysis to determine the equilibrium fluid-fluid demixing binodal, defined by equality of 
pressures and chemical potentials in coexisting colloid-rich and colloid-poor fluid phases. 
We restrict our study to parameters for which crystallization occurs only at higher colloid 
densities, well separated from fluid-fluid demixing. Figure El presents the resulting phase dia- 
grams for a fixed polymer-to-colloid size ratio, q = 1, and for effective polymer-to-colloid size 
ratios q' = 1, 0.8, 0.6, 0.4, corresponding to electrostatic colloid size ratios £ = 1, 1.11, 1.25, 
1.43, i.e., varying ionic strengths [varying k in Eq. (fTj)]. We assume, as discussed in Sec. IIII| 
that £ is independent of colloid density, which is valid in the high-salt-concentration regime 
(p s ^> Zp c ), where salt ions overwhelm counterions and k [Eq. Q] is essentially independent 
of p c . In the extreme limit, k — > oo (£ — > 1), we recover the binodal for the purely entropic 
case of neutral colloids. With decreasing ionic strength (decreasing k, increasing effective 
colloid radius), the binodal shifts to significantly higher polymer concentrations and the crit- 
ical point to slighty lower colloid concentrations. The qualitative consequence of lowering 
ionic strength is thus a significant enhancement of stability against fluid-fluid demixing. In 
the absence of polymer (j] p = 0), the mixture reduces to a one-component system of hard 
spheres of effective radius R' c . In this limit, crystallization occurs at effective colloid volume 
fraction r( c = £ 3 7/ c = 0.494, which is well above the demixing critical volume fractions for £ 
values considered here (see Fig. EJ). Note that the case £ = 1.43 is close to the limit beyond 
which demixing becomes metastable with respect to the freezing transition. At higher values 
of £, freezing may preempt demixing - a scenario that we do not explore here. 

The enhanced stability of charged-colloid - neutral-polymer mixtures can be physically 
interpreted in terms of a weakening of polymer depletion-induced colloidal attraction. Elec- 
trostatic repulsion increases the average separation between macroions, tending to lower the 
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frequency of configurations in which polymer is depleted from the intervening space. In our 
model, as the range of electrostatic repulsion increases (e.g., by removing salt), the effec- 
tive polymer-to-colloid size ratio decreases, reducing the range of the effective colloid-colloid 
attraction, and thus diminishing the driving force for phase separation. 

A physically equivalent interpretation of the electrostatic suppression of demixing follows 
from considering the polymer free- volume fraction a. From inspection of Eq. (|14jh it is not 
immediately obvious how a varies with electrostatic colloid size ratio £, given that, with 
increasing £, the quantities (1 — T]' c ), A, B, and C all decrease, while x increases. As Fig. |U 
illustrates, however, with increasing electrostatic repulsion between colloids, the free volume 
available to the polymer coils is monotonically reduced, even though the colloid-polymer 
interaction is unchanged. The reason for this reduction is that more strongly repelling 
colloids remain more widely separated, allowing less excluded volume to be hidden within 
overlapping polymer depletion shells surrounding the colloids. The lower free- volume fraction 
reduces the polymer entropy in the colloid-rich phase, thereby suppressing demixing. 

As a check on the consistency of the above interpretations, we consider the liquid-vapor 
interfacial tension 7 as a function of the colloid volume fraction difference Ai] c between 
the coexisting phases. To calculate 7, we exploit previous density-functional (DF) theory 
result a 19 i 62 for the AO model, which have been found to be in reasonable agreement with 
simulatio n 23 ! 24 . By appropriately rescaling the DF predictions of Ref. [l9| (An c — > Ai] c /^ 3 
and /37<t 2 — > /?7<7 2 /£ 2 ), we obtain the results shown in Fig. Evidently, the interfacial ten- 
sion increases with increasing effective colloid radius, which is consistent with the predicted 
suppression of liquid-vapor demixing. Although the mapping of DF theory from neutral 
to charged systems is straightforward in the absence of external potentials, more care may 
be required in the presence of a surface. For example, a hard wall for the neutral system 
becomes a non-additive wall for the charged system. In such cases, coupling to the external 
potential may have to be explicitly built into the theory. 

V. SUMMARY AND CONCLUSIONS 

In summary, we have investigated fluid-fluid demixing in mixtures of charged colloids and 
neutral non-adsorbing polymers. Our theoretical approach involves first mapping a model 
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mixture of charged hard-sphere colloids - interacting via an effective Yukawa electrostatic 
pair potential - and ideal polymer spheres onto an effective Asakura-Oosawa model with 
purely excluded-volume interactions, and then applying free-volume theory and density- 
functional theory to study, respectively, bulk phase separation and interfacial tension. The 
effective model is characterized by a single dimensionless parameter, namely the ratio of the 
polymer diameter to the effective colloid diameter. 

Within this simple framework, we find that increasing the range of the electrostatic 
repulsion between macroions (e.g., by decreasing salt concentration) suppresses demixing 
and increases the interfacial tension between coexisting colloid-rich and colloid-poor phases. 
Hence, our main conclusion is that electrostatic interactions stabilize the suspension against 
depletion-induced phase separation. As the depletion-induced attraction is short-ranged, the 
principal mechanism is analogous to charge-stabilization against coagulation induced by van 
der Waals forces. In both cases the electrostatic repulsion keeps the colloidal particles (at 
least partially) outside the range of attraction. Within our free- volume approach, the effect is 
clearly manifested through a reduction of the effective polymer size, and a correspondingly 
weaker depletion-induced attraction. The predicted increase of miscibility with decreas- 
ing salt concentration is qualitatively consistent with the observed behavior of mixtures of 
charged proteins and uncharged non-adsorbing polysaccharides (e.g., gelatin and dextran), 
which are constituents of many food colloids^. Future experiments and simulations would 
help to further test the qualitative trends predicted here. 

A possible extension of the present theory, which could yield more quantitative predic- 
tions, would be based on a variational method for the free energy. This method would split 
the free energy into two parts: a (zeroth order) reference term, describing the effective AO 
model, and a (first-order) perturbative term, describing the repulsive tail of the colloidal pair 
potential. Minimizing the total free energy with respect to the effective colloid hard-sphere 
diameter would give an upper bound on the true free energy. Such an approach might be 
better suited to describing colloid-polymer mixtures with long-range electrostatic repulsion 
between the colloids, which have been studied in recent experiment s^ 4 -^, and could be ap- 
plied also to mixtures of colloids and nonideal polymers, e.g., polyelectroryte a 40 i 41 i 42 i 43 i 44 , 
properly incorporating soft polymer-polymer interactions 29 . 

In this paper, we have focused attention on bulk fluid-fluid demixing in the high-salt- 
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concentration regime. As an outlook for future work, we mention two related phenomena that 
could be addressed within a similar theoretical framework. First, the fluid-solid (freezing) 
phase behavior of colloid-polymer mixtures is likely to be enriched by combining electrostatic 
and polymer-depletion interactions. It is well known^& that concentrated suspensions of 
charged colloids crystallize into close-packed fee or hep structures at higher ionic strengths, 
and into the more open bec structure at lower ionic strengths. It would be interesting to 
examine the competing influences of electrostatics and polymer depletion on freezing and 
the relative stabilities of crystal structures. 

Second, strongly deionized suspensions of highly charged colloids are known to exhibit 
surprisingly complex phase behavior. For example, experimental observation o 67 ' 68 and the- 
oretical prediction o 48 i 50 i 52 indicate the possibility of unusual bulk separation into macroion- 
rich and -poor phases, driven by a competition between counterion entropy and macroion- 
counterion attractive energy. Incorporating the volume energy [Eq. (jHJ)] into the total free 
energy would make possible an exploration of the influence of added depletants on counterion- 
facilitated phase separation in deionized suspensions. 
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FIG. 1: Sketch of the two models of charged colloids (c) and (effective spherical) neutral polymer 
(p). The solid curves indicate the actual sizes of the particles (radii R c and R p ) in the original 
model, and the dashed curves the rescaled sizes (radii R' c and R' p ) in the effective Asakura-Oosawa 
model. Note that the depletion layer thickness is the same in both models. 
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FIG. 2: Effective polymer-to-colloid size ratio q' = R' p /R' c vs. (a) actual polymer-to-colloid size 
ratio q = R p /R c and (b) electrostatic colloid size ratio £ = R' c /R c . 
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FIG. 3: Fluid-fluid demixing binodal of charged-colloid - neutral-polymer mixtures as a function of 
colloid and polymer volume fractions, rj c and rj p , for fixed polymer-to-colloid size ratio, q = R p /R c = 
1, and varying electrostatic colloid size ratio £ = R' c /R c . From bottom to top, £ = 1, 1.11, 1.25, 1.43, 
corresponding to increasing effective colloid radius. Symbols represent the respective critical points. 
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FIG. 4: Polymer free- volume fraction a [from Eq. (j!4j) ] as a function of colloid volume fraction r\ c 
for fixed polymer-to-colloid size ratio, q = R p /R c = 1, and varying electrostatic colloid size ratio 
£ = R' c /R c . From top to bottom, £ = 1, 1.11, 1.25, 1.43, corresponding to increasing effective colloid 
radius. 
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FIG. 5: Interfacial tension 7, in reduced units, between liquid (colloid-rich) and vapor (colloid- 
poor) phases vs. the difference in colloid volume fraction Ar] c between colloidal liquid and vapor 
(at bulk coexistence). 
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